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ABSTRACT 


The Marine Gas Turbine control systems in present use in the U. S. 
Navy are of such significant technological age that new design 
techniques could lead to more optimum performance and increased 
plant efficiency. To this end, a new real time Marine Gas Turbine 
simulation method is needed for advanced controller design and 
implementation. 

A modeling method is shown which utilizes real time sequential 
linearizations to approximate the true nonlinear response of the NPS 
Boeing 502-6A test facility. A validation of this simulation approach is 
presented. The method has immediate application to advanced 
controller design, especially to the design of modern regulators 
(Linear Quadratic), model reference controllers, and real time 
diagnostics. 
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I. INTRODUCTION 


Conventional control theory is limited to time invariant single-input- 
single-output systems and is usually based on a_ frequency domain 
design approach. Modern control theory, on the other hand, can be 
well applied to multiple-input-multiple-output, non-linear gas turbine 


systems which are now in use and it relies on a much more applicable 
time domain approach.! The feasibility of modern control theory 
such as Linear Quadratic Regulator (LQR) theory, in the application of 
marine gas turbines for propulsion has been demonstrated by past 
work at the Naval Postgraduate School in Monterey, California.2-4 For 
this application to be successful a simple, accurate, and robust real 
time marine gas turbine simulation must be developed. 

This thesis introduces a modeling approach to be applied in 
developing a real time simulation of engine response to meet the 
requirements of simplicity, accuracy, and robustness. It is simple in 
that it is only first order dependent, accurate in that it provides 
acceptable (+/- 10%) response for use in application of advanced 
engine control design, and robust in that it works well over the entire 


operational range of the gas turbine system. 


Il. PAST WORK IN GAS TURBINE MODELING 


Past work in Gas Turbine Modeling is discussed in terms of two 
areas: first, a brief discussion on gas turbine modeling in general is 
presented. This is followed by a discussion on the previous work 
completed at the Naval Postgraduate School in Monterey, California. 

Since most marine gas turbines have been derived from aviation 
sources, most of the early gas turbine modeling work was 
accomplished in the aviation arena. In the 1970's the Naval Gas 
Turbine Ship Propulsion Dynamics and Control Systems Research and 


Development Program was started to address the issues of designing 
Marine Gas Turbines for future use in the Navy.2 This included the 
dynamic modeling of the response of a marine gas turbine. Prior to 
this program, most of the work had been steady state modeling. 
However, the Navy realized that dynamic modeling was needed for 
advanced controller design. The program was successful in developing 
a gas turbine machinery dynamics and control system data base. From 
this work a computer-based propulsion control testbed was developed 
for use in gas turbine controller development. Cut and try linear 


control design then ensued, which resulted in a gain adjusted 


controller design in use in today's fleet. This approach resulted in a 
linear Proportional-Integral theory being applied to a nonlinear 
system. The results were deemed adequate for the Navy's uses. 

A new approach to this modeling is to develop a method to 
accurately model the propulsion plant nonlinear responses with a real 
time simulator. This would help eliminate or alleviate problems 
associated with the varying propeller loads felt in the gas turbine by 
predicting the responses before they actually occur, thus allowing 
more acurate control. More accurate control would mean better 
engine response to these loads and would, in turn, mean less wear, 
thus leading to longer engine life as well as better engine fuel 
efficiency. This is the direction taken at the Naval Postgraduate 
School, (NPS), in the gas turbine modeling problem. 

The previous work started with Phillip Johnson in 1985 in which 
he conducted the development and implementation of a load control 
system ( e.g. a water brake dynamometer) to emulate the shaft and 
propeller of a marine gas turbine.© The next step was taken by James 
Roger also in 1985. He used a simple multiport diagram and the 


Continous System Modeling Program (CSMP III) to incorporate an 


improved model for the dynamometer unloading valve into the existing 


system dynamic model.’ These two works formed the base for the 
next step performed by Vincent Herda in 1986. Herda developed a 


better multiport model and used it to develop a linear state space 


model and a_ nonlinear dynamic propulsion model.2 Robert Miller 
followed in the same year using Herda's models and examined the 
feasibility of using Linear Quadratic Regulator (LQR) control design 


theory for a small gas turbine, specifically the Boeing 502-6A test bed 


installation at NPS.3 Finally Vincent Stammetti in 1988 developed the 
first real time dynamic simulator for the gas turbine and compared an 


LQR controller to a classical Proportional Integral (PI) regulator based 


on his results with the simulator.4 To summarize, Herda validated 
cause and effect through the use of his new multiport model, shown in 
figure 1 , Miller validated the steady-state model method and 
Stammetti followed with the first linear model structure. It is 
important to note that these models were only accurate for one run 
of simulation and were not applicable over the whole range ol 
operation of the gas turbine. Their main point however, was to 
validate and show the feasibility of developing a better model for use in 
modern controller design to update those in current use today. To 
this end they were successful. This paper will continue this work and 


attempt to develop a simple, robust, and accurate real time marine gas 
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turbine simulator to be used for advanced controller design. 


il. MODEL STRUCTURE AND COARSE TUNING 


A. Modeling 

The gas turbine plant to be modeled was the 175 horsepower 
Boeing 502-6A which was assembled in a test bed configuration at the 
Naval Postgraduate School in Monterey, California. It was coupled to a 
Clayton 17-300 water brake dynamometer which was used to simulate 
the shaft and propeller loading on a marine gas turbine. The gas 
turbine/water brake assembly is shown in figure 2. The gas turbine 
itself was composed of two aerodynamically coupled sections, a gas 
generating section and a power output section. The gas generator had 
three major components. These were the compressor, the burner, 
and the high pressure turbine shown in figure 3. Connected by 
gearing off the gas generator shaft was the accessory section 
consisting of the fuel pump, lube oil pump, governor, tachometer 
generator, and the starter. The compressor was a single stage 
centrifugal compressor and was coupled to the high pressure turbine 
aerodymanically via two through-flow type, cross-connected 
combustion chambers. The high pressure turbine was a single stage 
axial flow high pressure turbine. It was aerodynamically connected to 
the low pressure or free power turbine which was also an axial flow 
turbine. The free power turbine output shaft was mechanically 


coupled to the water brake dynamometer which absorbs the energy of 
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Figure 2. Test Bed Configuration. 
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the free power turbine in simulating the action of a drive shaft and 
propeller. 

The modeling began with Herda's cause and effect multiport 
diagram shown in figure 1. The importance of this diagram to the 
method can hardly be overstated. Each of the components, the 
compressor, high pressure and free power turbines, and the gas 
generator and dynamometer shafts, were modeled by describing their 
outputs in terms of their inputs in equation form. For example, the 
compressor inputs are gas generator speed (Ng) and the high pressure 
turbine inlet pressure (P2) and the outputs are mass flow rate of air 
(Ma), high pressure turbine inlet temperature (T2), and compressor 
torque (Qc). The compressor block diagram is shown in figure 4. 

Ma 
Compressor T2 


P2 Figure 4. Compressor Block Diagram. 





Ng ¥Oc 


It is important to note that these input/output quantities as well as 
those for the other components must be measurable or calculable. 
They were recorded and or calculated over the operating range of the 
gas turbine (shown in figure 5), and the results are tabulated in table 


1. For clarity, the operating region of the turbine is shown in figure 5. 
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Figure 5. Gas Turbine Operating Range. 
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1.33353 
1.67667 
0.33933 
0.67333 
1.00667 
1.34967 


Ns 


$00 
1003 
1504 
2002 
2267 

$11 
1009 
1500 
2008 
2509 

S06 
1006 
PS@2 
2000 
2510 

S09 
1010 
1504 
2000 
Zo 

S09 
1010 
1510 
2017 
2510 


0.56818 
0.93182 
-31818 
» 18182 
.09091 
-81818 
65909 
~§2273 
- 36369 
. 18182 
. 11364 
-90909 
77273 
-59091 
-495455 
-45455 
pete? 5 
.04545 
-88636 
a hetion 
-90909 
-65909 
-40909 
- 25000 


=e eae Odrrr OO CO 0OrF0 0 Ee cdgqoca 


0.72371 
0.71561 
0.71573 
0.71573 
0.71573 
0.85760 
0.89969 
0.89975 
0.89975 
0.84148 
1.01243 
0.99619 
0.98411 
0.98411 
0.99825 
1.14700 
1.12666 
1.11886 
1.11484 
1.11993 
1.28258 
1.27065 
1.26275 
1.27544 


Ody iz F2 T4 P4 ME Ma QOc/ 
Ofpt Ohpt 
rpm - rpm-ftlbf£--f -- psig - f--psig-lbm/hr-lbm/s-ftlbf 

125 146.73 6.63 668 .85 7659857 s5cco 7 oce 

95 146.68 6.68 882 .94 78.39 1.71 60.981 

70 146.90 6.70 880 .96 79.40 1.71 62.764 

40 147.15 6.75 876 .95 79.40 1.71 63.326 

20 146.93 6.60 6878 .90 J9n40 517-7] 64 Cze 
180 177.45 9.18 877 1.15 "91045 2705 657628 
145 177.20 9.25 902 1.30 94°49 2.03 67) 4a 
115 176.15 9.40 902 1.35 94.98 2503569 222 

60 175.30 9.35 895 1.35 994798 2703 70.422 

40 174.98 9.33 897 1.30 94.98 2.01 72.434 
248 208.90 12.20 886 153390108. s1 2042757240. 2 
200 206.58 12.55 922 T9725 t1O Gs 2 foster 
170 208.03 12.43 917 1.80 )03 26052735074 
130 207.40 12.558 921.1. 90.13 Ceo oe 
100 208.50 12.55. 926 JSe5els. GU eZee oom 
320 244.18 16.05 913 2.10 123567 2.74 85,802 
280 247.33 16.05 960 2.20 131.44 2.69 80.095 
230 245.85 16.35 967 2.4907 135.976 2.67 978) 710 
195 244.36 16745 970 2.50 136n47e2-66 817167 
160 242.40 16.40 956 2.50 137. 24 2.66585. 736 
420 280.65 20.50 978 2. JON1S/28453-506 27772 .0 
36S 283.03 20.80 1002 mes T6i) 665s. Cs 94acr> 
310 281.05 20.90 1039 3.00 165.11°53502 903223 
275 263.50 20.95 1024 35308 le7760 33-04) 9 lmoan 
240 284.70 20.85 1044 3.30 168.34 2.99 89. 646 

(A) 

t2 p2 t4 p? mf ma qc/qhpt 
0.68247 0.947357 0.90795 0.449737 0.62841 0.73305 0.75705 
0.68223 0.97714 0.92259 0.499474 0.63992 0.72458 0.76226 
0.68326 0.97857 0.92050 0.50526 0.64816 0.72458 0.78455 
0.68492 0.498214 0.91632 0.50000 0.64816 0.72458 0.79157 
0.68390 0.97193 0.91841 0.47368 0.64816 0.72458 0.80025 
0.82535 0.65571 0.91736 0.60526 0.74653 0.86864 0.82120 
0.82419 0.66071 0.94351 0.68421 0.77135 0.86017 0.84319 
0.81930 0.67193 0.94351 0.71053 0.77535 0.86017 0.86661 
0.81535 0.66786 0.93619 0.71053 0.77535 0.86017 0.88049 
0.81386 0.66643 0.93828 0.68421 0.77535 0.85169 0.90542 
0.95767 0.87143 0.92678 0.83684 0.88416 1.02542 0.91035 
0.96084 0.89643 0.96444 0.90526 0.90473 1.00847 0.89222 
0.96758 0.88786 0.95920 0.94737 0.92738 0.99876 0.92797 
0.96465 0.89643 0.96339 1.00000 0.92735 0.99576 0.98382 
0.96977 0.89643 0.96862 0.97368 0.92735 1.00000 0.94584 
1.13572 1.14643 0.98502 1.10526 1.05037 1.16102 1.07252 
1.15037 1.19643 1.00418 1.15789 1.07298 1.13983 1.00119 
1.14349 1.16786 1.01151 1.26316 1.10588 1.13136 0.98389 
1.13656 1.17500 1.019464 1.31579 1.11909 1.12712 1.01459 
1.12744 1.17143 1.00000 1.31579 1.12033 1.12712 1.07170 
1.30835 1.46428 1.02301 1.42105 1.28849 1.29661 1.29762 
1.31642 1.48571 1.04812 1.50000 1.32131 1.28390 1.18356 
1.30721 1.49286 1.08682 1.57895 1.34784 1.27542 1.12776 
1.31860 1.499643 1.07113 1.73684 1.36816 1.28819 1.14542 
1.32419 1.98928 1.09205 1.73684 1.37420 1.26695 1.12057 


1.25077 


1.67333 


1.09091 


(A) Raw data 
(B) Normalized data 


1.25486 


Table 1. Steady State Data Taken Over Operating Range 
(Ng 20k -32k rpm, Ns 500 - 2500 rpm) 


ie 


At this point it should be noted that the data of table 1A was 
normalized according to each of the variable's midpoint for ease of 
understanding each variable's relative impact with respect to the 
others. Once normalized, the variables were are denoted using lower 


case letters. The resulting equation forms are: 


Ma = fi (Ng, P2), (1) 
T2 = fo (Ng, P2), and (2) 
Qc = f3 (Ng, P2). (3) 


The normalization values are shown in table 2. 


Table 2. Normalization Points. 


Ng = 26000 Ns = 1500 

Qd/Qfpt = 220 T2 = 215 
P2 = 14 T4 = 956 
P4= 1.9 Mf = 122.5 


Ma= 2.36 Qc/Qhpt = 80 


The next step was to distinguish between steady and dynamic 
components. Once this is complete, the component modeling can 
proceed. The rest of the component equation forms are shown in 


table 3. 
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Table 3. Component Equations. 


Gas Generator Shaft: Ng = fa (Qhpt, Qc) (4) 
Dynamometer Shaft: Ns = fs (Qfpt, Qd) (5) 
Free Power Turbine: P4 = fe (Ns, T4, Maf) (6) 
Qfpt = f7 (Ns, T4, Maf) (7) 

High Pressure Turbine: E = fg (Mf) (8) 
Maf = fo (Ma, Mf) (9) 

T4 = fio (Ma, T2, P4, Ng, E) (10) 

P2 = f11 (Ma, T2, P4, Ng, E) (11) 

Qhpt = fi2 (Ma, T2, P4, Ng, E) (12) 


B. State Selections 


State variables are the smallest set of variables which are required 
to ascertain the state or condition of a dynamic system.! The 
knowledge of these variables at an initial point in time (to), when 
combined with the input as time proceeds (t > to), will completly 
describe the systems behavior. The state of a system at any time (t) is 
indepentent of any prior state and input which precede the initial 
time point (to)!. Here, the turbine and load is the dynamic system and 


the states, when properly applied, will accurately describe the 
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dynamic behavior of the system. Significant prior research has 


determined that only three states are necessary to sufficiently 


describe the operating condition of the NPS gas turbine.® In reality 
however, only two states are necessary (ng, and ns) since the third 
state (e) represents a short first order lag of the input fuel (mf). This 
time lag is due to the burner action delay between the input of 
additional fuel and the thermodynamic energy output felt in the 
turbine. 

It is our goal to relate all of the dependent variables (p2, 
p4,t2,mf,...) to the three state variables ng, ns, and e to allow us to use 
the state space equation x = Ax + Bu (ref. 9) to describe the transient 
behavior of our system. Here x is the state vector and is composed of 
state variables ng, ns, and e, u is the input vector and is composed of 
mf and qd. Refering back to figure 1, note that we have redrawn 
Herda's system boundary at qd in order to have a more measurable 
input. Again, from prior research, it was known that each of the states 
is associated with a significant time lag in the system response. In 
this way, we knew that the fuel input component, the rotor shaft, and 
the output shaft must be regarded as dynamic components. We also 
knew that the remaining components could be regarded as 
instantaneous (memoryless), or steady components. The steady and 


dynamic components are shown in table 4. 
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Table 4. Steady and Dynamic Components. 


Steady Components: 
compressor ma t2 qc 
high pressure turbine p4_ aqfpt 
power turbine t4 p2 qhpt maf 


Dymamic Components: 


rotor shaft ng 
load shaft ns 
burner lag e 


C. Steady Components 

While the steady components of the model were given above as non- 
linear functions, in this section, we reduce these to linear 
expressions, where the coefficients were regressed and varried 
according to the engine operating point. 

The steady component equation form selection was accomplished 
through the use of an equation fitting program (Minitab), and plotting 
of each of the input variables versus an output variable minus the 
residual. If the residual plot displayed a parabolic shape, it was an 
indication that the input variable function in question was of a higher 
order. In the case of the equations with only two input variables a 


different method was used to fit the data. That is, one at a time, each 
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of the input variables were held constant and the other was plotted 
versus the output variable. Both methods are developed further in 
Appendix C. This resulted in accurate data fit equations. 

To simplify the equation reduction methodology, the steady 


component equations were kept in the form: 
x=cl y+c2z (13) 


For example, the compressor equations appeared as follows: 


t2=cl ng*+c2 ng+c3 p2+c4 (14) 
ma =c5 ng*+c6ng+c7p2+c8 and (15) 
qc = c9 ng? +clO ng+cll p2+cl2 (16) 


where the coefficients c] thru cl2 were the values from the equation 


fits. It is important to note the strong dependency of the state ng. 
This will appear again in the decoupling equations developed at the 
end of this chapter. The next phase of development was the 


formulation of the dynamic component models. 


D. Dynamic Components 

The dynamic components are derived from the gas generator shaft, 
the dynamometer shaft, and the lag relationship between mf and e. 
The dynamic components were modeled with differential equations 


which, when integrated over time will determine the state of the 


Lalé 


system. The dynamic components took the form: 


ng = qhpt-qe (7) 
ns = qfpt-qd (18) 
and e = (mf-e)/T. (19) 


These, along with the previous equations of the steady components, 
were next linearized to form a dynamic equation set for the simulation 


model. 


E. Linearized Dynamic Equation Set 

The linearized dynamic equation set was developed from the 
perturbational or derivative values of the steady state equations 
(symbolically, the perturbational variables will be proceeded by a "d"). 
The compressor equations are used below to demonstrate the 
formulation of the dynamic equations. Recall that the t2 steady state 


equation was: 


t2 =cl ng? + c2 ng+c3 p2+c4 (14) 
So, the dynamic equation took the form: 
dt2 = [ dt2/dng ] dng +[dt2/dp2]dp2 (20) 
Here the partial differentials will be denoted with upper case letters 


so that the resulting form of equation 20 is: 


dt2 = Al dng+A2 dp2. (21) 


The rest of the steady and dynamic equations previously developed 
were also converted to dynamic form in this manner. The results 


appear below in table 5. 


Table 5. Dynamic Equation Set. 


Compressor at2="Al-dng + A2 dp2) San. (22) 
dma = B] dng + B2 dp? ............ (23) 
dqe—"@ wdne+4C2 dp2 ........... (24) 
Gas Generator Shaft 
JGB dng = dqhpt-dac ..............cceeceee (25) 
Dynamometer Shaft 
JOB dns.— dotot-dad avmenee ts. (26) 
Free Power Turbine 
dp4 = D1 dns+D2 dt4+D3 dmaf ....... (27) 
dqfpt = F1 dns+F2 dt4+F3 dmaf ..... (28) 
High Pressure Turbine 
Gea (CimimMecle) / i) wuemameine, ... :..Snences cs (29) 
from continuity-- dmaf = K1 dmf + K2 dma................ (30) 
dt4 = G] dma + G2 dt2 + G3 dp4 + G4 dng + G5 de ............ (31) 
dp2 = H1 dma + H2 dt2 + H3 dp4 + H4 dng + H5 de ........... (32) 
dqhpt =11 dma+ 12 dt2 +13 dp4 +1!4dng +15 de .......... (33) 


These resulting equations from table 5 will be used to develop the 
state space relationships through a reduction methodology to get 


three equations ( one each for ng, ns, and e) in terms of ng, ns, e, mf, 
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and qd. 


F. Decoupling the Dynamic Equation Set 

In the past, during the equation set reduction, a number of the table 
3. equations were used more than once. This resulted in the 
development of singularity points in the simulations. To avoid this, 
another method was developed in the present work. This was 
achieved by inspecting the system block diagram (fig. 1) and 
deciding that if two equations, one for p2 and the other for p4, 
could be expressed in terms of the states ng, and ns only, then it 
would be possible to reduce the equation set directly (without 
resubstitution) thus eliminating the singularities. The new p2 and p4 
equations were developed with the assistance of the Minitab 
regression routine and steady state data for ng and ns. This resulted 
in the conclusion that both variables had a second power effect with 
ng and a first power effect with ns. With these relationships in mind, 


the computer routine was used to fit the data in a two step process. 


First , equations for p2 and p4 were generated in terms of ng2, ng, 
and ns. Next, the data generated from these equations for p2 and p4 
was plotted versus ng at five different constant ns values. These plots 
were compared to the actual data to check for accuracy of the 
equations. The results are shown on the two graphs in figure 6 and 


show the robustness of these equations over the entire range of the 
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Figure 6. P2 and P4 graphical Representations. 
Note: Solid lines are equations, symbols are data. 
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data. The new dynamic equations for p2 and p4 are now in the form of: 
p2=cl ng42+c2ng+c3ns+c4 and (34) 


p4 =c5 ng*2 + c6 ng + c7 ns + C8. (35) 
These equations were then linearized into the form: 


dp2 = Hi] dng + H2 dns foreqn. 11 and (36) 


dp4 = Dj dng + D2 dns for eqn. 6. (37) 
These twelve equations were then reduced with the reduction 
methodology as described in detail in Appendices A and B. 

The results of the reduction transformed the linearized equations to 
state space equations, which were used in the simulation model. The 
resulting state space equations from Appendix B are: 

dng = Z1 dng + Z2 dns + Z3 de 
JGB JGB JGB 


dns = Z4 dng + Z5 dns + Z6 de + Z7 dmf - 1 dqd 
JDB JDB JDB JDB JDB 


de = - 1de+_1dmf 
T T 
where Z1 through Z7 were derived from combinations of the Al 


through K2 coefficients of the dynamic equations. The standard state 


space equation is x = A(t) x +B(t) u.9 Here x is a matrix made up of 
dng, dns, and de, x is a state vector composed of the state variables: 


dng, dns, and de, and u is a vector composed of inputs dmf and dqd. 


mee 


The A and B matrices are: 


all = 2.1. Al2Qi=22 WN 5) = 7A! 
JGB JGB JGB 
A21 = Z4 A22 =Z9 A23 = Z6 
JDB JDB JDB 
A31 =O A32 =O A33 =-L 
T 
and Bl1=O0O B12 =O 


EB? l= 27 B22 =-l. 


JDB JDB 
B31 = _L B32 =O 
ae 


and the state space equations now appear as: 


dng = All dng + Al2 dns + Al3 de + B11 dmf + B12 dmf 


dns = A2] dng + A22 dns + A23 de + B21] dmf + B22 dmf 


de = A31 dng + A32 dns + A33 de + B31 dmf + B32 dmf. 


A coarse tuning was next performed to determine the best 
numerical values of the coefficients. This was necessary since the first 
set of regressed coefficients did not provide accurate simulation. This 
was attributed to the coarse instrumentation that is present on the 
NPS turbine and the manner of fitting in the regression routine. 
Consequently, attention was focused on the following coefficients as in 
need of tuning: Al, A2, Bl, B2, Cl, C2, Fl, andI1 -15. Those for Dl, 


D2, Hl, H2, Kl, and K2 were found previously using the data 
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regression program. Initially chosen for a starting point was that for 
B2 (dma/dp2). This was an obvious choice for much is known about 
the relationship between ma and p2 in compressors. To find B2, p2 
was plotted versus ma at constant ng values. From this the slope 
oma/dp2 (B2) was obtained. The Bl expression (@ma/dng) was found 
by regression using the plotted value of B2 in the Minitab routine. A2 
and C2 were then found using plotted values to find the slopes dt2/dp2 
(A2) and dqc/dp2 (C2). These were then used as inputs into the 
relationships: 


Bl=dma/dng = Al=dt2/dng = Cl =dqc/dng 
B2 dma/dp2 A2 dt2/dp2 C2 dqc/dp2 


from which Al and Cl were obtained. The results were then checked 
and verified through the use of the data regression program. At steady 
state qc = qhpt is required, thus the relationship for Cl (dqc/dng) = 14 
(dqghpt/dng). So, with 14 determined, the value of I] (dqhpt/dma) can 
be found as I1 = 14/B1 and [2 (dqhpt/dt2) is 14/Al. With these three, 
I1, 12, and 14, as inputs, I3 (dqghpt/dp4) and I5 (dqhpt/dmf) were 
obtained through the use of the data regression routine. The 
numerical values can be found in Appendix C. 

The next step involves fine tuning the computer simulation model 
and results in the setting of the remaining coefficients of the A and B 


matrices. 
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IV. FINE TUNING THE MODEL 


In fine tuning, the actual dynamic data was used in conjunction with 
a simulation approach to select final values for the model constants. 

The simulation approach used the simple linear state space 
equations for dng, dns, and de. These described a first order 
differential equation set, with time varying A and B matrices, to mimic 
the actual nonlinearity of the system response. This was accomplished 
by approximating the true nonlinear relationships with a series of 
small linear steps as depicted in figure 7. 

The dynamic simulation computer model was developed from a 
routine written by Herda and later modified by Stammetti, it is shown 
in figure 8. The major change to the previous program was in the A 
and B matrix equations. In the past, the equations were complex 
and were a combination of exponentials and second order equations. 
In the present model they were simple, first order linear functions of 
the states. Another change which occurred in the fine tuning process 
was the adjusting of the dynomometer shaft inertia. This change 
was for shaping the response and did not affect the steady state 


response of the model. 


A. Computer Simulation 


The simulation program was built using four sections. These were 
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State 








Sequential a 


Linearizations 






Actual Response 


Time 


Figure 7. Sequential Linearization. 


Input: U 
A(t), B(t) 


Compute: x = A(t) x + B(t) u 


X= f x dt 


Figure 8. Simulation Process. 
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the dynamic data input, the initial condition, the dynamic, and the 
derivative sections. The actual dynamic data was taken on a four 
channel strip chart recorder and depicts the real turbine response of 
Ng, Ns, Qd, and Mf as a function of time. These dynamic runs were 
compiled in Appendix E and provided the source and data for the 
simulation. The initial condition section established the initial 
conditions for each particular run. The dynamic section formulated 
the display of the dynamic data and was followed by the derivative 
section which formulated the simulated response. The outputs of 
these final two sections were then plotted for comparison. These 
plotted results were used to assist in a grid search on the still doubtful 
component coefficients to determine the final values of the A and B 
matrices. The procedure was to guess values of the coefficients, then 
simulate the model response. If the model and the data agreed, then 


the search was over. 


B. The Grid Search 

A grid search technique was applied in two and three parameter 
searches to find the optimal values for the remaining coefficients. It 
was formed using a the spread sheet formatted program called Excell, 
by Microsoft software. An example of this is found in Appendix C. The 
top portion provided the inputs for the coefficients used in the 


calculation of the A and B matrices. The middle portion calculated the 
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intermediate Y variables. These were then combined with the 
previously mentioned A thru K coefficients to calculate the A and B 
matrix values, three of which were ng dependent. 

We also knew that the steady state required x = 0. This fact was 
used to generate additional equations which must be satisfied by any 
coefficient set to guide the grid search. The last section of the spread 
sheet provided for the calculation of the values of dng and dns for the 
simulation runs 3 and 9 as well as the peak dns value for run three 
(figure 9). Run three, shown in figure 9, developed into the most 
difficult run for determining the best coefficients for Ns simulation 
and is the reason for the addition of the peak dns value. The values 
for the G, F2, and FS3 coefficients were the focus of the grid search as 
we had least confidence in them. These are the important coefficients 
in determining the value of dns and in turn, the ns simulation values. 
The G1 - GS values were calculated in the same manner as described 
eariler for the I1- I5 values. They represent the partial derivatives of 
the equation for t4 (eqn. 10) and were determined to have little effect 
on the A21, A22 and A33 results. Their elimination as a search 
parameter reduced the search to a two parameter one. Various values 
for F2 and F3 were tried, and the resulting A and B matrix coefficients 
calculated were inserted in the computer simulation. The results 
were plotted on the NPS mainframe computer system, using the 


easyplot routine, to check for accuracy of simulation. This process was 
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Figure 9. Dynamic Data Runs. 


repeated until the optimal F2 and F3 values were obtained. 
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V. VALIDATION OF THE MODEL 


The validation of the model was a two step process which carried 
over from the fine tuning. After the fine tuning was complete, and 
the model resulted in the proper shape for runs 3 and 9, and with the 
steady state response within the required tolerances, the first step 
was complete. The next step was to take the model and apply it 
elsewhere in the operating range of the turbine to check for accuracy. 
If the proper response in shape and prescribed tolerance criteria (+/ 
- 10 percent in the dynamic and steady response) were still met, 
then the validation was completed. This process was only performed 
around the borders of the operating range (figure 9). If successful, 
then the model should be applicable anywhere in the turbine 
operating range. 

The first run modeled was the one for run nine which models the 
Ns transient response from 500 to 2500 rpms while holding the value 
for Ng, as load varies, constant at 32000 rpms (+/- 5%) for the 
conditions shown in run nine in Appendix E.. Figure 10 is the 
graphical representation of the modeled run and the dynamic data run 
for this transient. The criteria has been met for run 9, so the model 
was applied to the validation run one. Run one, simulates the 
transient for Ns from 500 to 2215 rpms (the turbine will not operate 


above 2215 rpm at low Ng) with the value for Ng constant at 20000 
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Figure 10. Simulation-Run Nine. 
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Figure 11. Validation-Run One. 
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rpms (+/- 5%); the results are shown in figure 11. Again, the 
accuracy criteria were met. This covers the upper and lower bounds 
of the operating range. 

Next, run three was used in conjunction with run nine in the fine 
tuning portion. This run is a vertical run on the left side of the 
operating range and is an Ng transient from 20000 to 32000 rpms 
with Ns being held constant at 500 rpms to conform with the actual 
dynamic data run. This run is shown in figure 12. Note that the 
accuracy criteria are easily met in Ng but apparently not satisfied in 
Ns. However, when the accuracy criteria is applied to the normalized 
Ns, the transient is acceptable. To validate this run the model is next 
applied to the right hand side of the operating range. This is run 
seven and, again, is a Ng transient from 20000 to 32000 rpms with Ns 
slightly increasing from the lower constraint of 2215 rpm to 2500 
rpms which is attainable at the upper limits of the operating range. 
This run is presented in figure 13. A validation run was attempted 
vertically in the center of the operating range. This is run five and 
it's response was of the same result shown in figure 11 for run seven 


and therefore was not shown as a separate figure. 
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Figure 12. Simulation-Run Three. 
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Figure 13. Validation-Run Seven. 
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wt NCLUSION 


A modeling method has been developed for real time simulation of a 
marine gas turbine. The results of this modeling method should be 
used with the LQR controller design theory to develop a new 
controller for the Boeing test bed at NPS. This controller should then 
be tested to ascertain it's ability to be more fuel efficient than the 
present system. 

This model was simple in that it used matrix elements for A and B 
which were only linear state functions. The dependency was only on 
the Ng state. This simplicity will provide fast run times on a small 
computer for model reference control and diagnostics. 

Accurate simulation performance has been demonstrated in the 
compressor/high pressure turbine section and the resulting steady 
states have been achieved to +/- 3%. An exception to this accuracy is 
the constant Ns simulation for the vertical runs . This is only an 
indication that the final values for F2, F3, and possibly G1 are not fine 
tuned enough to provide for more accurate constant Ns response. 
Investigation was conducted into this irregularity to determine the 
cause of the steady state results being to high. The fact that the steady 
state response for a constant Ns is higher in runs five and seven was 
probed and a number of solutions were looked at. One such solution 


was to check for the possibility that the oscillating torque for runs five 
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and seven shown in Appendix E was causing the Ns value to continue 
to rise and settle out at a higher steady state. These oscillations do 
not appear to be the cause nor have any affect on the steady state 
response of Ns. This was verified by modifying the computer 


simulation routine so as to remove the oscillation. Another solution 
was to adjust the A2] values until a proper steady state value for ns was 
reached. At the time of publication this has not provided any 
satisfactory results. Another possibility is that the A21 coefficients 


may have to change the sign to produce the proper response. The last 
area to examine as a possible solution to this problem is to input a 
varying inertia for the water brake dynamometer. In the past this 
inertia has been assumed to be a constant value when, in fact, it 
actually changes with the addition and subtraction of water. 

The last recommendations for this project encompass the entire 
project. First, a digital data acquisition system should be used to get 
the steady state and dynamic data. This will make the fine tuning 
process with the data regression program more accurate. Finally, one 
computer system should be used for the whole process, from the 
regression and grid search to the computer simulation (vice the two 
system used in this project). In this way, results could be achieved 
and correlated in a more efficient manner. 

The model was robust in that, with the exception of the constant 


Ns in the vertical runs, the model was applicable to all areas of the 
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operating range. This was supported by the excellent response in Ng 


throughout the whole operating range. This shows that the 
coefficients for the All, Al2, and Al3 were accurate and the 
procedure to find them was sound. 

Once the proper values for A21, A22, and A23 are found and this 


modeling technique has been successfully demonstrated on the Boeing 
engine this method will have real time application to the future 
controller design for the General Electric LM- 2500 and the follow on 


marine gas turbine controllers designed for the future. 
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APPENDIX A 


COMPONENT MODELING EQUATIONS : 


Compressor 
dt2 =A] dng+A2ap2 ........ yi 
dma = B] dng + B2dpz2 _.............. 2 
dqc =C] dng+ C2dp2 “227 3 


Gas Generator Shaft 


qhpt JGB dng = dghpt-dcc ...... 


Dynamometer Shaft 


afpt qd JDB dns = dafpt-dad ...... 
ns ns 


Free Power Turbine 


maf Goa =D) dng + D2 dns ............ oi 

t4 

p4 dqfpt = F1 dns+F2 dt4+F3 dmaf... 7 
qfpt } ns 


High Pressure Turbine 


de =(dmf-de)/T _...... 8 
mda maf 
t2 t4 
p2 p4 





From Continuity: 


dmaf = K] dmf + K2 dma....... 9 


dt4 = G} dma + G2 dt2 + G3 dp4+G4dng+G5de_................... 10 
Bret Mtl CLG GU cape) cess crerenaesteccecvesssssccncseatecsavsscnsenertesenceeess 1 
dghpt =1]1 dma+ 12 dt2 +13 dp4 +14 dng +1[5de ................... 12 


Note: 1. Lower case variables represent normalized variables. 


2. The d _ represents the perturbation of the variable. 


* Decoupling equations. 
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APPENDIX B 
EQUATION SET REDUCTION METHODOLOGY: 


The equations in appendix A are reduced into the state space 


components of the A and B matrices. 


Substitute eqn 11: 


dp2 = H] dng + H2 dns 


et CME aclein vans daiiosevcesaesvavsececeveds 1] 
into eqns 1, 2, and 3 for dp2. 
eee OP ine + (ADHD) GS o..c.s.csiccecee.ccesccecccsscoscescoosceecesscosees l 
ieee ae lets 1) Ci + (2D) ANS ..........ccccceeceosscossssnesccovsscvcccvsessees 2 
Salome nde Dil) ING © (COHOAMS <..5....ccsssscceesssccccesescoccovscccceveveoreeee 3 
In eqn 12: 
dqhpt = 1] dma + 12 dt2 +13 dp4+1l4dng+iI5de ............... 12 
substitute eqns 1, 2, and 6 for dt2, dma, and dp4. 
dghpt = 11 (B] + B2H1) dng + 11 (B2H2) dns + 
I2 (Al + A2H]1) dng + I2 (A2H2) dns + 
13 (D1) dng + 13 (D2) dns + 
14 dng Se ls oss ctle be av dete soca 12 
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Collect terms: 


dqhpt = (11 (B] + B2H1) + I2 (Al + A2H1) + [3 (D1) + 14) dng + 


(11 (B2H2) + I2 (A2H2) + 13 (D2) ) dns + I5 de 


Let Yl = I1 (B1 + B2H1) + I2 (Al + A2H]1) + 13 (D1) + 14 
Y2 = 1] (B2H2) + I2 (A2H2) + 13 (D2) 


Then: dqhpt=Yldng+Y2dns+I5de .........eee 12 


Take eqn 4 and substitute in new eqn 12 for dqhpt and eqn 3 for 
dqc. 


JGB dng = Y1 dng + Y2 dns + I5 de - 


(C1 + C2H1) dng - (C2H2) dns 
Collect terms: 


dng = ( (Yl - (C1 + C2H1)) dng + (Y2 - (C2H2)) dns + I5 de )/JGB 


Let: . Zl = (Ve ieee o a) 
ZZ = 12 - (C22) 
Z3 = 15 

Then: 


dng = (Z1 dng +Z2dns+Z3de)/JGB  —_____n eee seen Rl 
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In eqn 10: 
dt4 = G] dma + G2 dt2 + G3 dp4+G4dng+G5de _............. 10 
substitute eqns 1, 2, and 6 for dt2, dma, and dp4. 
dt4 = G1] (B] + B2H]1) dng + G1 (B2H2) dns + 
G2 (Al + A2H1) dng + G2 (A2H2) dns + 
G3 (D1) dng + G3 (D2) dns + 


G4 dng IGE CME rE reer 10 
Collect terms: 


dt4 = (G1 (B1 + B2H1) + G2 (Al + A2H1) + G3 (D1) + G4) dng + 


(G1 (B2H2) + G2 (A2H2) + G3 (D2) ) dns + G5 de 


Let Y3 = G1 (Bl + B2H]1) + G2 (Al + A2H]1) + G3 (D1) + G4 
Y4 = G1] (B2H2) + G2 (A2H2) + G3 (D2) 


Mem, G64 =—YSdng+ Y4dnSt+ G5 dE ————aecesecaceccccececscesess 10 


In the continuity eqn (9) substitute eqn 2 for dma. 
dmaf = K] dmf + K2(B1 + B2H1) dng + K2(B2H2) dns _.......... 9 
Substitute new eqn 10 for dt4 and new eqn 9 for dmaf into: 
eqn 7: daqfpt = Fl dns + F2 dt4 + F3 dmaf 


This results in: 
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dqfpt = F1 dns + F2Y3 dng + F2Y4 dns + F2G5 de + 
F3K1 dmf + F3K2(B1 + B2H1) dng + F3K2(B2H2) dns 


Collect terms: 
daqfpt = (F2Y3 + F3K2(B1 + B2H1)) dng + 


(Fl + F3K2(B2H2) + F2Y4) dns + F2G5 de + F3K1 dmf 


Let: Z4 = F2Y3 + F3K2(B1 + B2H1) 


Z9 = Fl + F3K2(B2H2) + F2Y4 


Z6 = F2G5 
Za el 
Then: 
dqfpt = Z4 dng +Z5dns+Z6de+Z7dmf un... sane 7 


Substitute new eqn 7 for gfpt in eqn 5. 
JDBdns=dqfpt-dqd eee eee ee eeeeee eee 5 


dns = (24 dng + Z5 dns + Z6 de + Z7 dmf - dqd)/JDB _................. R2 


Finally: Eqn 8 is R38. Rl, R2, and R3 make up the state space 


equation set. 


46 


The state space equation set is the following: 


dng=Zl dng+Z2dns+Z3de —_—sC_—_. nurssssssssveeseeees Rl 
JGB JGB JGB 
dns = Z4 dng + Z5 dns + Z6 de + Z7 dmf- 1 ddqd _....... R2 
JDB JDB JDB JDB JDB 
de = - J}de+ Jdmf na... eee. eee R3 
T T 


The resulting A and B matrices are: 


Alge= Zh Ao = Z2 Als = 23 
JGB JGB JGB 


A2] = Z4 A22 = Z5 A23 = Z6 


JDB JDB JDB 
A31 =O A32 =O A33 = a) 
il 


and 
Bi l= © B12 =O 


B21 =Z7 B22 =-1 


JDB JDB 
B31 = _1L B32 =O 
T 
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APPENDIX C 


PARAMETER SELECTION METHODOLOGY: 


The parameter selections were accomplished through three means: 
1. Plotting actual slopes from data. This was shown 
in chapter III section F. 
2. Data reduction using mainframe program 
Mini-Tab. 
3. Grid searching using coefficients from 1. and 2. 
and known relationships which must be satisfied. 
The spread sheet is from the Micro-soft program 
Excell and was used on an Apple Macintosh 
computer. 

The data reduction program Minitab was used for all of the 
regressions performed. The input for the regression program was the 
steady state turbine data shown in table 1 A. This routine was used to 
determine the order of the independent variables in the equations for 
Ma, T2, and Qc, as well as to regress the expressions for P2, P4, T4, 
Qhpt and Maf. To determine the order of a independent variable a 
three step process was followed. First, the dependent variable was 
regressed using the data for the independent variables. In doing this a 
residual was calculated for each set of the independent variables. This 


was then subtracted from the regressed value of the dependent 
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—— = 


variable. The resulting difference was plotted versus each of the 
independent variables. In the plot for the regression of T2, shown in 
figure C-la, there is a definite parabolic shape which is an indication 


that the independent variable in question (Ng) is of a higher order. To 


verify this, T2 was regressed again, this time using Ng42, Ng and P2. 
The plotting was repeated with the resultant plot shown in figure C- 
lb. This time the parabolic shape was not in evidence and the 
predictor percentage, which indicates the accuracy of the regression 
based on the data, increased from 99.6 to 99.9 percent indicating an 
excellent regression for the t2 equation. The resulting partials for all 
of the component expressions are shown in figure C-2 lines 3 through 
21. 

The spreadsheet grid, figure C-2, was used in conjunction with the 
simulation program to search out the best values for the rest of the 
coefficients in the equations in appendix A. These were input for Gl, 
F2, and F3 in lines 3-21 and then used in lines 25-32 to calculate the 
values of the coefficients of the A and B matrices. Lines 34-38 show 
the proper signs for these matrix coefficients as well as the dns and 
dng relationships which must be satisfied during the search process. 
Once these were Satisfied, plots were generated and checked for the 


accuracy of the simulation. 
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Figure C-1. Regression Plots 
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Figure C-2. Sample Grid Search Spread Sheet. 
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APPENDIX D 


COMPUTER SIMULATION CODES: 


The computer codes were written using the Dynamic Simulation 
Language (DSL). 


Simulation Run 1: DSL code for dynamic Ns run at constant Ng 
(20000 rpm). 


Simulation Run 3: DSL code for dynamic Ng run at constant Ns 


(500 rpm). 


Simulation Run 5: DSL code for dynamic Ng run at constant Ns 
(1500 rpm). 


Simulation Run 7: DSL code for dynamic Ng run at Ns 


( 2252-2500 rpm). 


Simulation Run 9: DSL code for dynamic Ns run at constant Ng 


(32000 rpm). 


Sf 


HHH HK KKK KKK KKK KKK KKK KK KKK KKEKKEKEKKEKEKEKKKKKKKKKKKKKKKKKEKKKKKKKKKKEK 


* 

* 

BOEING MODEL 502-6A GAS TURBINE * 

* 

DYNAMIC COMPUTER SIMULATION * 

* 

MODIFIED BY * 

S.D. METZ * 

FROM ROUTINES BY * 

V.J.HERDA AND V.A. STAMMETTI * 

* 

THIS PROGRAM SIMULATES THE DYNAMIC RESPONSE OF THE NPS * 
BOEING GAS TURBINE TEST FACILITY USING A MULTILPLE * 
LINEARIZATION TECHNIQUE. * 
* 
RAKKKKKEKKKEKKEKEKEKKEKEKEKEKKKKEKKKEKKEKEKKEREKKEKKKEKKEKKKKKKKKKKKKKEKKEKKKKKKKKKKKKKkKK 


ARAM JG=0.009525, JD=0.3000, PI=3.14159, TC=1.0 


THE FOLLOWING VALUES LISTED ON THE FUNCTION CARD ARE FOR FUEL FLOW, 
GAS GENERATOR SPEED, TORQUE AND DYNO SPEED AS A FUNCTION OF TIME. 
THESE VALUES WERE OBTAINED FROM STRIP CHART RECORDS AND ARE ENTERED 
Pr tHE FORM (E.G. FUEL FLOW) ...TIME( SEC), FUEL FLOW..... 


THIS SET IS FOR EXPERIMENTAL RUN # 1. 


QANQNAA * * * e+ ee £ION © & FH He HH He He + HE He ee 


BGeN NGDATA= 0.0,20230.0, 5.0,20230.0, 10.0,20230.0, 15.0,20230.0, 
Bene coz Ono, 2o.0),20230.0, 50.0,20230.0 

MeemhetioDALA>= ©.0,505.0, 1.0,505.0, 2.0,536.7, 3.0,568.33, .. 
ae OO Ome O oS). 7,.0.0,695.0, 7.0,790.0, 8.0,885.0, 
eo Oss eee tO.1O, 1170.0, 11.0,1328.3, 12.0,1486.7, ... 
omOpee OCs mel asO,1eo5. 0, 15.0,1930.0, 16.0,1993.3, 
iWeOrZOso wy elo. 0 2050.3, 19.0,2120:0, 20.0,2151.7, 
POMZeOr omeee. O17) .0, 23,0,2183.3, 24,0,2215.0, 
wo Os. OO. 072215.0 

SeeeeltirDATA= 0.0,78.39, 5.0,78.39, 10.0,78.39, 12.0,79. 89, 
oO om. orm 7.0,00.2060, 190,781.38, 20.0,82.50, 
WeeOnoc elope eon ,82.13, 50.0,82.13 


eee OODATA= ©.0,125.0, 1.0,125.0, 2.0,118.3, 3.0,115.0, 
Aeon ogy, 105.0, 6. 0,98.33, 7.0,85.0, . 
One Tey meso, oe. 33, 10.0,45.0, 11.0,31. 67, 
PaO eeOg elias, 1.670, 13.0,5.0, 15.0,5.0, 
Zon; 5.07 50. 0,5. 0 


* 
* 
* 


INITIAL 
“8 ESTABLISH INITIAL CONDITIONS. 
* 


NGI=20230. 0 
NSI=505.0 
MFI = 78.39 
QDI = 125.0 

* 

* 

* 

* SET INITIAL STATE PERTURBATION TO ZERO 
Pnca= 0. 0 
DNS = 0.0 
DEE) 0. 0 
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NGB 
NSB 
MFB 


QDB 


OCB 
JGB 


JDB 
DYNAMIC 
* 


* DAT 


* 


NGD 
NSD 
MED 
QDD 
MFDLY 
bs STAT 
All = 
Al2 = 
Al3 = 
A21 = 
A22 = 
A23 = 


A31 


il 


A32 


A33 = 


Bi 


B12 
B21 °=. 
B22 = 
B31 = 


B32 = 
* 
a 


DERIVATIVE 
NOSORT 
* 


* COMPUTE 

* RUN #9 
DMF 
DOD 


DNGDOT 
DNSDOT 
DEDOT 


+ + 


26000. 0 
1500.0 

122.5 

220.0 

80.0 
(NGB/QCB)*JG 
(NSB/QDB) * JD 


A CURVE FORMULATION 


AFGEN( NGDATA, TIME) 
AFGEN( NSDATA, TIME) 
AFGEN( MEDATA, TIME) 
AFGEN( QDDATA, TIME) 
TRANSP( 100,MFI,0. 70,MED) 


E SPACE LINEAR MODEL FORMULATION 
(20. 282-32. 560*( NGL/NGB) )/( JGB) 
(-0. 4480) /( JGB) 

(13. 624)/( JGB) 

(2. 1131+. 91786*( NGL/NGB) )/( JDB) 
(-. 0998-0. 445*( NGL/NGB) )/( JDB) 
(-0. 1885) /( JDB) 

0.0 

0.0 

-1.0/TC 

0.0 

0.0 

035250/( JDB) 

-1.0/( JDB) 

1207506 


O.0 


INPUT TO THE NONLINEAR MODEL, ME(T), WW(T). 


=(MEDLY - MFI)/MEFB 

=(QDD = QDI) /QDB 

Al1*DNG + A12*DNS + A13*DE + B1l1*DME 

A21*DNG + A22*DNS + A23*DE + B21*DMF + B22*DQD 
A33*DE + B31*DME 


* DYNAMIC EQUATIONS FOR LINEAR MODEL. 
DNG=INTGRL( 0. 0, DNGDOT) 
DNS=INTGRL( 0. 0, DNSDOT) 

DE =INTGRL( 0. 0,DEDOT) 


NGL 
NSL 
SORT 


= NGI + DNG*NGB 
= NSI + DNS*NSB 
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THE STATEMENTS IN THE PREVIOUS (DERIVATIVE) SECTION YIELD VALUES 
OF 'NG', AND 'NS' AS CALCULATED BY THE NONLINEAR AND STATE SPACE 


MODELS. THE STATEMENTS BELOW COMPUTE THE VALUES OF 'NG', AND 'NS' 
AS RECORDED FROM GAS TURBINE TEST DATA. 


+e t+ + 


* 


TERMINAL 

METHOD RKS 

CONTROL FINTIM=50. 0, DELT=0. 001 
PRINT 0.5,NGD,NGL,NSD,NSL,MED,QDD 
END 

S.OP 
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(Re RRR ERE RRR RRR ER RRERRRR RRR RRR RRR RRR RRR RRR RRR RRR ERR RRR ERP CeCe 


* 

* 

BOEING MODEL 502-6A GAS TURBINE * 

* 

DYNAMIC COMPUTER SIMULATION s 

* 

MODIFIED BY * 

S.D. METZ * 

FROM ROUTINES BY * 

V. J. HERDA AND V.A. STAMMETTI 

* 

THIS PROGRAM SIMULATES THE DYNAMIC RESPONSE OF THE NPS * 
BOEING GAS TURBINE TEST FACILITY USING A MULTILPLE * 
LINEARIZATION TECHNIQUE. 
* 
KEKE KKKKKEKEEKKAEEKEEKEKEEKEEEEEKEEEEKEZEEREEZKEEKEREZELEZELEEEKEEEEEKEKEXKEKERS 


ARAM JG=0.009525, JD=0. 3000, PI=3.14159, TC=1. 00 


THE FOLLOWING VALUES LISTED ON THE FUNCTION CARD ARE FOR FUEL FLOW, 
GAS GENERATOR SPEED, TORQUE AND DYNO SPEED AS A FUNCTION OF TIME. 
THESE VALUES WERE OBTAINED FROM STRIP CHART RECORDS AND ARE ENTERED 
IN THE FORM (E.G. FUEL FLOW)  ...TIME( SEG ye eUE cee ee 


x 
. 
. 
* 
® 
& 
« 
* 
® 
« 
. 
« 
* 
. 
* 
* 
Cc 
Cc 
r 
+ 
4 
. 
« 
: 
. 
¢ 
Cc 
Cc THIS SET IS FOR EXPERIMENTAL RUN #% 3. 
C 
Cc 
AEFGEN NGDATA= 0.0,20220.0, 3.55,20220.0, 3.57 32030527 732 c7 20 2 ear 
4.0,20698.4, $.0,23186.1, &. 0725482547 O72 eee 
8.0,28926.9, 9.0,30649.1, 10.0, 3208050 Os 2257 ee 
11.0,32180.0, 11.5,32084.3, 12.0) 32h etye eo eee 
AFGEN NSDATA= 0.0,506.0, 2.0,505.7, 3.6,/S05" Sete oe o 
5.0,508.3, 6.0,509.7, 7.0,511. 1,78) 0pe ee ee 
10.0,514. 4, 10.4,514.7, 11.0,/514°54 Siva, 
13.0,515.8, 14.0,516.7, 16:3,517227 eee eee, 
18.0,515.8, 19.0,514.9, 20,060,513. 0, Sea eee ee 
AFGEN MEDATA= 0.0,77.5, 3.0,77.5, 353,79. 6)c pee oer 
4.0,91.99, 5.0,106.5, 6.0,118.9, 720 ie Ges co. 0o) 2 eee 
9.0,156.2, 10.0,172.7, 11.0,156. 2, V2Z5Oyteer ee eee 


AFGEN QDDATA= 0.0,125.0, 3.6,125.0, 3.8,123.4,5537070)) 
5.0,172. 4, 6.0,226.7, 7.0,274%. 1, S207 32s seer ocr 
10.0,430.0, 10.5,443.6, 11.0, 4300, 12407412 
13.0,416.4, 14.0,419.8, 15.0,416.4,) te5@ 4222307 
20.0, 430, 0 (50. 0, 43¢ 


* 
* 


INITIAL 
* ESTABLISH INITIAL CONDITIONS. 
* = 
NGI=20220. 0 
NSI=506. 0 
MEI = 77.5 
ODile=) 1250 
* 
* 
* 
* SET INITIAL STATE PERTURBATION TO ZERO 
* 
DG =i0n0 
DNS = 0.0 
Deo .0 
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NGB = 26000.0 
NSB = 1500.0 
MFB = 122.5 
QODB = 220.0 
QCB = 80.0 
JGB = (NGB/QCB)*JG 
JDB = (NSB/QDB)*JD 
DYNAMIC 
* 
* DATA CURVE FORMULATION 
* 
NGD = AFGEN( NGDATA, TIME) 
NSD = AFGEN( NSDATA, TIME) 
MED = AFGEN( MEDATA, TIME) 
QDD = AFGEN( QDDATA, TIME) 
MFDLY = TRANSP( 100,MFI,0. 70,MED) 
* 
* STATE SPACE LINEAR MODEL FORMULATION 
* 
All = (20. 282-32. 560*(NGL/NGB) )/( JGB) 
* 
Al2 = (-0. 448)/(JGB) 
* 
A13 = (13.624)/(JGB) 
* 
A21 = (2.1131+.91786*(NGL/NGB) )/( JDB) 
* 
A22 = (-. 0998-0. 445*(NGL/NGB) ) /( JDB) 


A23 =(-.18851)/(JDB) 


A31 = 0.0 
* 

A32 = 0.0 
* 

Reg=)=1,0/7TC 
* 

Bll = 0.0 
* 

B12 = 0.0 


B21 =.035250/( JDB) 
B22 = -1.0/(JDB) 
B31 = 1.0/TC 


B32 = 0.0 


* 
* 


DERIVATIVE 
NOSORT 
* 


* COMPUTE INPUT TO THE NONLINEAR MODEL, ME(T), WW(T). 
* RUN #3 


DMF =(MEDLY - MFI)/MEB 
DQD =(QDD - QDI)/ODB 


DNGDOT = All*DNG + A12*DNS + A13*DE + B11*DMF 
DNSDOT = A21*DNG + A22*DNS + A23*DE + B21*DME + B22, 000 
DEDOT = A33*DE + B31*DME 


* 


* DYNAMIC EQUATIONS FOR LINEAR MODEL. 
DNG=INTGRL( 0. 0, DNGDOT) 
DNS=INTGRL( 0. 0, DNSDOT) 

DE =INTGRL( 0. 0,DEDOT) 
NGL = NGI + DNG*NGB 
NSL = NSI + DNS*NSB 


SORT 
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m THE STATEMENTS IN THE PREVIOUS (DERIVATIVE) SECTION YIELD VALUES 

bi OF 'NG', AND 'NS' AS CALCULATED BY THE NONLINEAR AND STATE SPACE 

si MODELS. THE STATEMENTS BELOW COMPUTE THE VALUES OF 'NG', AND 'NS' 
= AS RECORDED FROM GAS TURBINE TEST DATA. 

& 

TERMINAL 

METHOD RKS 


CONTROL FINTIM=50. 0, DELT=0. 001 

PRINT 0.50,NGD,NGL,NSD,NSL,MED,QDD, DNSDOT,A21,DNG, 
A22,DNS,A23,DE,B21,DMF,B22,DQD 

END 

STOP 
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* * 
* * 
* BOEING MODEL 502-6A GAS TURBINE * 
* * 
* DYNAMIC COMPUTER SIMULATION * 
* * 
MODIFIED BY * 
* S.D. METZ * 
* FROM ROUTINES BY * 
* V.J.HERDA AND V.A. STAMMETTI * 
* 
THIS PROGRAM SIMULATES THE DYNAMIC RESPONSE OF THE NPS * 
BOEING GAS TURBINE TEST FACILITY USING A MULTILPLE * 
LINEARIZATION TECHNIQUE. 
* 
kKkk&k&kkk&kkkkkkk&kkkkkkikkkkkkkkkkkkkkkkkkkkkkkkdkkkkkkkkkkkkkkkkkkkkkkkkk 


ARAM JG=0.009525, JD=0.3000, PI=3.14159, TC=1.00 


* 

* 

* 

* 

* 

* 

e 

c 

P 

& 

* THE FOLLOWING VALUES LISTED ON THE FUNCTION CARD ARE FOR FUEL FLOW, 
* GAS GENERATOR SPEED, TORQUE AND DYNO SPEED AS A FUNCTION OF TIME. 
* THESE VALUES WERE OBTAINED FROM STRIP CHART RECORDS AND ARE ENTERED 
* IN THE FORM (E.G. FUEL FLOW) ...TIME(SEC), FUEL FLOW..... 

* 
* 
E 
Cc 
Cc 
g 
A 


THIS SET IS FOR EXPERIMENTAL RUN # 5. 


BeemeNGDATA> O.0,20220.0, 3.1,20220.0, 3.6,20412.1, 4%.0,21372.6, 
Seen zsb99.7,) 6.0,25086.9, 7.0,27711.8, 8.0,29440.6, . 
PmemoOl77 ao. O,32206.8, 10.2,32322.1, 10.8,32130.0, 

ioe sO. Oem i2.1,3920148.7, 13.4,32206.8, 14.5,32014.7, 
Pemoroc2oCne,n lo. 6,32014.7, 17.0,32206.8, 50.0,32130..0 

* 

AFGEN NSDATA= 0.0,1491.6, 1.0,1490.4, 3.0,1489.6, 4.0,1492.5, 
Soo Ol Ofcom 1506.5, 7.0, 1511.7, 7.8,1513.3, . 
Pao le oO, 1513.8, 12.0,1508.3, 13.0, 1510.0, 
5070, 1510. 0 

* 

BeeeveMEDATA= 0.0,83.75, 3.2,83.75, 5.0,109.6, 6.0,122.2, 
pre ees. oo. 0, 150.2, 9.0,162.8, 9.8,172. 4, 
MemOnwtoo.?, 1. 0,165.0, 50.0,165.0 

* 

AFGEN QDDATA= 0.0,71.0, 4.4,71.0, 5.0,78.7, 6.0,94.1, 
MeomlOvo, 6.0,420.9, 9.0,134.4, 10.2,149.8, ... 
omorian. 9, le 6,447.9, 12,4,143.2, 13.0,145.9, 
13.5,144.0, 14.0,145.9, 15.0,144.0, 50.0,144.0 


* 
* 
* 


INITIAL 


‘ ESTABLISH INITIAL CONDITIONS. 


*eFe & © 


NG1=20220. 0 


NSI=1491. 6 
MET = 83.75 
QDI = 71.0 
SET INITIAL STATE PERTURBATION TO ZERO 
DNG = 0.0 
DNS = 0.0 
DE = 0.0 
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NGB = 26000.0 
NSB = 1500.0 
MFB = 122.5 
QDB = 220.0 
OCB = 80.0 
JGB = (NGB/QCB)*JG 
JDB = (NSB/QDB)*JD 
DYNAMIC 
* 
* DATA CURVE FORMULATION 
* 
NGD = AFGEN(NGDATA, TIME) 
NSD = AFGEN( NSDATA, TIME) 
MED = AFGEN(MEDATA, TIME) 
QDD = AFGEN( QDDATA, TIME) 
MEDLY = TRANSP(100,MFI,0.70,MED) 
& 
* STATE SPACE LINEAR MODEL FORMULATION 


All = (20, 282-32. 560*( NGL/NGB) )/( JGB) 
Al2 = (-0. 4480) /( JGB) 

Al3 = (13.624)/( JGB) 

A21 = (2.1131+.91786*( NGL/NGB) )/( JDB) 
A22 = (-. 0998-0. 445*( NGL/NGB) )/( JDB) 
A23 =(-0.18851)/( JDB) 

NL = 1), © 

R32 = Oud 


AGS. = Sino 


Bll O20 


Biz |= s0"e 
B21 =.035250/( JDB) 
B22 = -1.0/(JDB) 
Boi =2 oy ae 


B32 = 070 
* 


* 


DERIVATIVE 
NOSORT 


* 
* COMPUTE INPUT TO THE NONLINEAR MODEL, ME(T), WW(T). 
* RUN #9 

DMF =(MFDLY - MFI)/MEB 

DOD =(QDD - QDI)/QDB 


DNGDOT = Al1*DNG + Al2*DNS + AlS*DE + Baa fie 
DNSDOT = A21*DNG + A22*DNS + A23*DE + BZi*DMES + 352270 
DEDOT = A33*DE + B31*DME 


+ 


* DYNAMIC EQUATIONS FOR LINEAR MODEL. 
DNG=INTGRL( 0. 0, DNGDOT) 
DNS=INTGRL(0. 0, DNSDOT) 

DE =INTGRL(0.0,DEDOT) 
NGL = NGI + DNG*NGB 
NSL = NSI + DNS*NSB 


CORT 
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* THE STATEMENTS IN THE PREVIOUS (DERIVATIVE) SECTION YIELD VALUES 
* OF 'NG', AND 'NS' AS CALCULATED BY THE NONLINEAR AND STATE SPACE 
* MODELS. THE STATEMENTS BELOW COMPUTE THE VALUES OF 'NG', AND 'NS' 
* AS RECORDED FROM GAS TURBINE TEST DATA. 

* 

TERMINAL 

METHOD RKS5 


CONTROL FINTIM=50. 0,DELT=0. 001 

PRINT 0.5,NGD,NGL,NSD,NSL,MFD,QDD,DNSDOT,A21,DNG,A22,DNS, 
A23,DE,B21,DMF,B22,DQD 

END 

STOP 
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* 
* 
BOEING MODEL 502-6A GAS TURBINE x 
* 
DYNAMIC COMPUTER SIMULATION * 
* 
MODIFIED BY * 
S.D. METZ * 
FROM ROUTINES BY * 
V. J. HERDA AND V.A. STAMMETTI * 
* 
THIS PROGRAM SIMULATES THE DYNAMIC RESPONSE OF THE NPS * 
BOEING GAS TURBINE TEST FACILITY USING A MULTILPLE * 
LINEARIZATION TECHNIQUE. * 
* 
KEKE KKKEKKKKKKKEKKEKKEKAKKKKKKKKKKEKEKKEKKEKKAKKEKKEKEEKKKEKEKKEKKKKKKEKSE 


ARAM JG=0. 009525, JD=0. 3000, PI=3.14159, TC=1. 00 


THE FOLLOWING VALUES LISTED ON THE FUNCTION CARD ARE FOR FUEL FLOW, 
GAS GENERATOR SPEED, TORQUE AND DYNO SPEED AS A FUNCTION OF TIME. 
THESE VALUES WERE OBTAINED FROM STRIP CHART RECORDS AND ARE ENTERED 
IN THE FORM (E.G. FUEL FLOW) ...TIME(SEC), FUEL FLOW..... 


* 
* 
* 
* 
« 
n 
* 
* 
~ 
* 
- 
* 
* 
* 
* 
* 
Cc 
Cc 
5 
* 
* 
* 
* 
* 
* 
* 
c 
Cc THIS SET IS FOR EXPERIMENTAL RUN # 7. 
C 
Cc 
AFGEN NGDATA= 0.0,20220.0, 4.0,20220.0, 4.5,2075927,, 5.0 (213527 
6.0,24460.3, 7.0,26541.9, 8.0,28122.4, 9.0, 23/760, ee 
10.0,31399.0, 10.4,32092.9, 10,673 2367 8 ye ee 
12.0,32170. 0, 15.0,32170.0;7 20507 25,0 

AFGEN NSDATA= 0.0,2252.0, 1.0,2245.4, 2 O00 220 Serre 20 eee 

4.0,2225. 6, 5.0,2278.4, 6.0,2436. 6373 2occee, 

7.7,2621.6, 8.0,2611.0, 9.0,23735. 47 eee oe 

11.0,2542.4, 13.0,2516.0, 15750) 25167Cye 2 
AFGEN MEDATA= 0.0,78.8, 1.0,78.8, 3.0, 78.68, 4.0757 ee 

4.5,95.3, $.0,103.6, 6.0,311 9) 7. OZ er oe oor 

9.0,161.7, 10.0,174. 1, 10.4,17853,)) iO 

15. O77 OOF 20.407 1 On 


AEFGEN QDDATA= 0.0,25.0, 2.0,25.0, 4.0, 25.073 520. 5n- 
6.0,49.6, 7.0,92.6, 8.0,129. 4, 970) ocr 
10.0,233.9, 11.0,264.6, 12.0, 25273 32 en- eee 
15. 0,,240. 0, 20707 22080 


* 
* 
* 


INITIAL 
* 


* ESTABLISH INITIAL CONDITIONS. 
NGI=20220.0 


Nol=2252-0 
MFI = 78.75 
ODI = 25.0 
* 
* 
* 
‘i SET INITIAL STATE PERTURBATION TO ZERO 
DNG = 0.0 
DNS = 0.0 
DE = 0.0 
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NGB = 26000. 0 
NSB = 1500.0 
MFB = 122.5 
QDB = 220.0 
QCB = 80.0 
JGB = (NGB/QCB)*JG 
JDB = (NSB/QDB)*JD 
DYNAMIC 
* 
* DATA CURVE FORMULATION 
* 
NGD = AEGEN(NGDATA, TIME) 
NSD = AFGEN(NSDATA, TIME) 
MED = AFGEN(MEDATA, TIME) 
QDD = AFGEN( QDDATA, TIME) 
MEDLY = TRANSP(100,MFI,0.70,MED) 
* 
* STATE SPACE LINEAR MODEL FORMULATION 


All = (20. 282-32. 560*( NGL/NGB) )/( JGB) 
Al2 = (-0. 448)/( JGB) 
Al3 = (13.624)/(JGB) 
A21 = (2.1131+.91786*(NGL/NGB) )/( JDB) 


A22 = (-. 0998-0. 445*( NGL/NGB) )/( JDB) 


A23 = (-.18851)/( JDB) 
a 

A31 = 0.0 
* 

A32 = 0.0 
* 

A33 = ~1.0/TC 
* 

B11 = 0.0 
® 

B12 = 0.0 


B21 =.035250/( JDB) 


B22 = -1.0/( JDB) 
* 
B31 = 1.0/TC 
* 
B32 = 0.0 
* 
* 
DERIVATIVE 
NOSORT 


* 


x COMPUTE INPUT TO THE NONLINEAR MODEL, ME(T), WW(T). 
™ RUN #7 
* 


DMF =(MEDLY - MFI) /MEB 
DOD =(QDD - QDI) /QDB 


DNGDOT = All*DNG + Al2*DNS + A1l3*DE + B11*DME 
Pew tee = n2zi*DNG + AZ2*DNS + A23*DE + B21*DMEF + B22*DQD 
DEDOT = A33*DE + B31*DME 


* 


* DYNAMIC EQUATIONS FOR LINEAR MODEL. 
DNG=INTGRL( 0. 0, DNGDOT) 
DNS=INTGRL( 0. 0, DNSDOT) 

DE =INTGRL( 0. 0,DEDOT) 
NGL = NGI + DNG*NGB 
NSL = NSI + DNS*NSB 
SORT 
* 
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* OF 'NG', AND 'NS' AS CALCULATED BY THE NONLINEAR AND STATE SPACE 
* MODELS. THE STATEMENTS BELOW COMPUTE THE VALUES OF 'NG', AND 'NS' 
* AS RECORDED FROM GAS TURBINE TEST DATA. 

* 

TERMINAL 


METHOD RK5 

CONTROL FINTIM=20. 0,DELT=0. 001 
PRINT 0.25,NGD,NGL,NSD,NSL,MED,QDD 
END 

STOP 
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* 

* 

BOEING MODEL 502-6A GAS TURBINE * 

* 

DYNAMIC COMPUTER SIMULATION * 

* 

MODIFIED BY * 

S.D. METZ * 

FROM ROUTINES BY . 

V.J.HERDA AND V.A. STAMMETTI * 

& 

THIS PROGRAM SIMULATES THE DYNAMIC RESPONSE OF THE NPS * 
BOEING GAS TURBINE TEST FACILITY USING A MULTILPLE * 
LINEARIZATION TECHNIQUE. * 
& 
KAKKKKEKKKEKKKEKKEKEEKKKKAKKKEKEKEKEEKEEKEEKEKKEEEKEKEKEKEKEEEKEKEKKEKEKEK 


ARAM JG=0.009525, JD=0.3000, PI=3.14159, TC=1. 00 


THE FOLLOWING VALUES LISTED ON THE FUNCTION CARD ARE FOR FUEL FLOW, 
GAS GENERATOR SPEED, TORQUE AND DYNO SPEED AS A FUNCTION OF TIME. 
THESE VALUES WERE OBTAINED FROM STRIP CHART RECORDS AND ARE ENTERED 
Mueine FORM (E.G. FUEL ELOW) ...TIME( SEC), FUEL FLOW..... 


THIS SET IS FOR EXPERIMENTAL RUN # 9. 


QAAAD * * &€ + ee HeYVOHO © + eee eee ee He He EH 


PeeemenolAtA= ©.0,32150.0, 5.0,32150.0, 10.0,32150.0, 15.0,32150.0, 
POmope 2) 2e1Opm25.0,32190.0, 50.0,32150.0 

Peete SDATA= 0.0,510.0, 2.0,510.0, 3.0,576.1, 4.0,642.14, ... 
Deon OSes, o.0,774.3, 7.0,840.4, 8.0,922.9, 9.0,1005.5, 
omnia etl, 1269.8, 12.0,1435.0, 13.0,1583.6, . 
oe oe oO, 1903.5, 16,0,2128.7, 17.0,2227.8, 
ieee nceZzcws eto. e, 2069.0, 20.0,2409.5, 21.0,2459.1, 
MemOmra oe lee On2o0G.6, 24.0,2525.0, 25.0,2525.0, 
SOM ,Zo25. 0 

eee eDATA= 0.0,155.8, 3.0,155.8, 5.0,155.8, 7.0,160. 0" 
eo ow, Ole. O, 164.1, 15.0,164.1, 18.0,164.1, 4... 
emo mOe no eco. O7os.3, 24.0,168.3, 25.0,168.3, 50.0,168.3 


Roem OUDATA= 0.0,394.0, 2.0,394.0, 3.0,382.5, 4.0,371.0, 
Porc o ee nooo. 5, §.0,325.0, 9.0,313.50, ... 
Poms oe sO met, 204.8, 12.0,261.8, 13.0,250. 3, 
iS mOmeoo MO moro, e2l.5, 16,0,221:5, 19.0,221.5, ... 
Pomeneenn eee ine 210.0, 923.0,210.0, 25.0,210.0, 50.0,210.0 


* 
* 
* 


INITIAL 
* 


- ESTABLISH INITIAL CONDITIONS. 
* 


NGI=32150.0 


NSI=510.0 
MFI = 155.8 
QDI = 394.0 
* 
* 
* 
* SET INITIAL STATE PERTURBATION TO ZERO 
* 
DIG = 0, 0 
DNS = 0.0 
Dew = 0.0 
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NGB = 26000.0 
NSB = 1500.0 
MFB = 122.5 
QDB = 220.0 
OCB = 80.0 
JGB = (NGB/QCB)*JG 
JDB = (NSB/QDB)*JD 
DYNAMIC 
* 
DATA CURVE FORMULATION 
& 
NGD = AFGEN( NGDATA, TIME) 
NSD = AEFGEN( NSDATA, TIME) 
MED = AFGEN( MEDATA, TIME) 
ODD = AFGEN( QDDATA, TIME) 
MEDLY = TRANSP(100,MFI,0. 70,MED) 
* 
* STATE SPACE LINEAR MODEL FORMULATION 


All = (20. 282-32. 560*( NGL/NGB) ) /( JGB) 


Al2 = (-0. 4480) /( JGB) 
*& 
Al3 = (13.624)/(JGB) 
* 
A21 = (2.1131+.91786*( NGL/NGB) )/( JDB) 


A22 = (-. 0998-0. 445*( NGL/NGB) ) /( JDB) 


A23 =(-0.18851)/( JDB) 


A31 = 0.0 
* 
Aa2 = 080 
*& 
A33 = -1.0/TC 
*& 
Bll = 0.0 
* 
B12 = 0.0 


B21 =.035250/( JDB) 
B22 = =10 0/0) 


B31 = 1.0/TC 


B32 0.0 


® 
* 


DERIVATIVE 

NOSORT 

* 

* COMPUTE INPUT TO THE NONLINEAR MODEL, ME(T), WW(T). 
* RUN #9 


DMF =(MEDLY - MFI1)/MEB 
DQD =(QDD - QDI1)/ODB 


DNGDOT = All*DNG + Al2*DNS + Al3*DE tee re 
DNSDOT = A21*DNG + A22*DNS + A23*DE + B2Z1*DME + B227 bap 
DEDOT = A33*DE + B31*DMF 


+ 


DYNAMIC EQUATIONS FOR LINEAR MODEL. 
DNG=INTGRL( 0. 0, DNGDOT ) 
DNS=INTGRL( 0. 0,DNSDOT) 

DE =INTGRL( ©. 0, DEDOT) 
NGL NGI + DNG*NGB 
NSL NSil + DNS“ NSB 


SORT 
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THE STATEMENTS IN THE PREVIOUS (DERIVATIVE) SECTION YIELD VALUES 
OF 'NG', AND 'NS' AS CALCULATED BY THE NONLINEAR AND STATE SPACE 
MODELS. THE STATEMENTS BELOW COMPUTE THE VALUES OF 'NG', AND 'NS' 
AS RECORDED FROM GAS TURBINE TEST DATA. 


+ + + + 


* 


TERMINAL 

METHOD RKS 

CONTROL FINTIM=50. 0,DELT=0. 001 

PRINT 0.5,NGD,NGL,NSD,NSL,MFD, QDD,DNSDOT,A21,DNG, 
A22,DNS,A23,DE,B21,DMF,B22,DQD 

END 

STOP 
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APPENDIX E 
SOURCE DATA RUNS: 


Run 1: Data run for dynamic Ns run at constant Ng. 


Ng -- 20230 - 20240 rpm Ns -- 505 - 2215 rpm 
Qd -- 125 - 25 Ibft Mf -- 78.4 - 82.1 lbm/hr 


Run 3: Data run for dynamic Ng run at constant Ns. 


Ng -- 20220 - 32180 rpm Ns -- 506 - 513 rpm 
Qd -- 125 - 480 Ibft Mf -- 77.5 - 160.3 lbm/hr 


Run 5: Data run for dynamic Ng run at constant Ns 


Ng -- 20220 - 32130 rpm Ns -- 1492 - 1510 rpm 
Qd -- 71 - 144 lbft Mf -- 83.8 - 165 Ibm/hr 


Run 7: Data run for dynamic Ng run at Ns (2252-2516 rpm). 


Ng -- 20220 - 32170 rpm Ns -- 2252 - 2516 rpm 
Qd -- 25 - 240 Ibft Mf -- 78.8 - 170 Ibm/hr 


Run 9: Data run for dynamic Ns run at constant Ng. 


Ng -- 32150 rpm Ns -- 510 - 2525 rpm 
Qd -- 394 - 210 Ibft Mf -- 155.8 - 168.3 lbm/hr 
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